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In this work we will advance farther along a line previously developed 

concerning our proposal of a time interval operator, on finite dimensional 

spaces. The time interval operator is Hermitian, and its eigenvalues are time 

values with a precise and interesting role on the dynamics. With the help of 

the Discrete Phase Space Formalism (DPSF) previously developed, we show 

that the time interval operator is the complementary pair of the Hamilto- 

nian. From that, a simple system is proposed as a quantum clock. The only 

restriction is that our results do not apply to all possible Hamiltonians. 



I. INTRODUCTION 

In quantum mechanics the concept of time has been an issue always submitted to con- 
tinuous scrutiny and discussion, specially if one aims to treat it as an observable. Although 
very interesting results have been already achieved, the question is still far from to be closed. 
As it was shown long ago by Pauli [1], direct attempts to find an operator which obeys a 
canonical commutation relation with the Hamiltonian and, at the same time, that have a 
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reasonable physical interpretation are bounded to fail. In fact, there are in the literature 
not too many proposals of Hermitian operators which can be somehow related to time. 

In this connection there are interesting approaches to the question of time as an observ- 
able in quantum mechanics which should be mentioned. For example, the method which 
takes benefit of a doubled Hilbert space [2-5], while others explore the possibility of defining 
observables as positive-operator- valued (POV) measures [6]. In what refers to the intriguing 
concept of time in the context of barrier penetration there are several different approaches, 
although most of them do not concern a time operator [7-9], exception made to [10]. 

Notwithstanding all this interesting approaches, there seems not to be an attempt which 
can, at the same time, display a known algebra satisfied by the Hamiltonian and the 'time 
operator', with a physical interpretation for all quantities involved, and with the time op- 
erator being an orthodox Hermitian observable. From another point of view, we must also 
mention the work by Pegg [11], which looks for a complementary operator to the Hamilto- 
nian, that have some results similar to those of this paper. In the following we shall develop 
an idea which we recently proposed [12], that fulfills all the conditions mentioned above 
provided that the physical systems under consideration satisfy given conditions. 

The starting point of our approach is based on the consideration of quantum degrees 
of freedom characterized by finite dimensional state spaces. As has been shown, in these 
cases it is possible to make use of the Discrete Phase Space Formalism (DPSF) as a suitable 
auxiliary tool [13-15] since this formalism can be seen as a generalization of the Weyl-Wigner 
formalism to degrees of freedom without classical counterpart. As such, operators acting 
on the appropriate Hilbert (ket) spaces can be mapped onto functions of integer indices 
provided that a basis in the operator space is well defined. By its turn, it has been proposed 
that the functional form of the operator basis elements consist, in general terms, of a double 
discrete Fourier transform of the Schwinger operator basis [16]. Once the action of the 
operator basis elements on relevant state kets is defined, the mapping procedure, that takes 
operators to the corresponding representative function in the discrete phase space, can be 
put to practical terms. In section II we briefly review the main ideas referring to the operator 
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basis. 

Since we dispose of a mapping procedure obtained from the DPSF, operator equations 
can be converted into ordinary equations involving the phase space counterparts (the mapped 
expressions) of the respective operators. If we restrict ourselves to degrees of freedom of the 
type we are considering, with no classical counterpart, this convenient feature permits us to 
show that it is possible to find an operator, which we call time interval operator, that is a 
generator of energy shifts and its exponential obeys a Weyl-Schwinger commutation relation 
with the time evolution operator. This particular operator, defined for such systems, is 
a Hermitian operator such that its spectrum is a discrete set of time intervals with an 
important role in the dynamical evolution. This is, however, only accomplished by a certain 
class of Hamiltonians, and the Weyl-Schwinger algebra is fulfilled only for some prescribed 
time intervals, which are in fact the eigenvalues of the time interval operator. In this 
connection we shall see that there is no general time interval operator but rather a different 
time interval operator for each particular Hamiltonian which admits such a 'canonical' pair. 
The time evolution governed by this particular time interval operator, satisfying that algebra, 
characterizes a particular dynamics since it evolves the state of the system over the sites of 
the discrete phase space occupying only one site (exhausting the total probability) at each 
time interval, and must not be confused with the general time evolution governed by the 
usual continuous time, which, by its turn, can always be carried out. This discussion is 
basically what is accomplished in sections III and IV. 

Although the proposal of a discrete time has been discussed in other contexts, in this 
work we are not making any hypothesis about the background c-number valued parameter t, 
but rather discussing the properties of an operator and its eigenstates, a set of discrete time 
intervals which we denote 'quantum clock times', as we propose a simple quantum clock in 
the end of section IV. The move from the independent variable time to the concept of time 
interval is the key to the whole process, that in Section V is presented with some general 
remarks and conclusions. 
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II. THE OPERATOR BASIS 



Considering a time independent Hamiltonian and its set of eigenstates consisting of an 
N— dimensional ket space, 



H\u k ) = E k \u k ), k = 0,1, ...N -1, 



(1) 



we construct the operator 



JV-1 

V = \ U k)( u k+l\, 
k=0 



(2) 



where hereafter we adopt a cyclic notation 



\Uk) = |«fc(modA0) 



(3) 



which in particular means that in Eq.(2) the last term of the summation is in fact \un-i)(uo\. 
From this definition it follows that 



V s \u n ) = \u n - s ). 



(4) 



Now the operator V can in turn be diagonalized,and once its eigenstatesare obtained, it 
can be directly verified that the well known Schwinger unitary operators algebra results are 
reproduced [16], that is, 



U\u k ) = exp 



2ni 
~N 



k 



K), V\v k )=exp 



2ni 



k 



\v k ), 



(5) 



u N = v N = i, 



(6) 



V S \u n ) = \Un- s ), U S \v n ) = \v n+s ), 



(7) 



being that in our case the set {|w„)} diagonalizes both the Hamiltonian and the Schwinger 
operator U. Furthermore, the two sets of states are connected by a discrete Fourier transform 



(v k \u n ) 



exp 



2ni 



N 



-kn 



(8) 



and the operators obey the Weyl-Schwinger algebra 

-2m 



U J V l = exp 



N 



■Jl 



(9) 



It has been also shown that the set of N 2 operators [15] 



JV-l 
2 



G(m,n) = — U j V l exp 



cxp 



2m 



(mj + nl) 



(10) 



constitutes a complete and orthogonal basis in operator space; the relevant properties of 
this basis are collected in reference [15], from which we recall the main results. In the form 
that it is presented, this basis is suited to deal with odd dimensional spaces, although it can 
be easily adapted to even ones. In fact, as this basis keeps, from the Schwinger basis, the 
property of decomposition into independent sub-bases, each with a prime dimension [16], 
we, for simplicity, shall hereafter consider only prime iV's, although most of the results in 
what follows do not depend on this condition. We are also omitting the so called modular 
phase, which should, in principle at least, be present in the basis, but is irrelevant in the 
particular present calculations. 

Thus, any operator acting on those finite dimensional state spaces can be written in the 
form 



iV-l 
2 



= — J2 o(m,n)G(m,n), 



(11) 



Ji 2 



where the coefficients (the operator representative on the discrete phase space, or equiv- 
alently, the discrete Weyl- Wigner mapped function associated to the operator) are given 
by 



o(m,n) = Tr[G\m, n)0\. 



(12) 



III. TIME EVOLUTION AS A LATTICE SHIFT 

Systems with a time independent Hamiltonian H have a time evolution operator of the 
form 
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G(At) = exp 



-iHAt 



h 



(13) 



that propagates the states \ip(t)), to which corresponds the density operator P(t) = 
\ip(t)) (ip(t)\ for a general pure state (although this is not an essential assumption for what 
follows) in such a way that 



P(t + At) =g(At)P(t)G\At) 



(14) 



Now let us map this operator equation onto the discrete phase space assuming the hy- 
pothesis 



G{Ar) = U-\ ke[0,N-l], 



(15) 



that is, the unitary time evolution operator is equal to some integer power of the Schwinger 
unitary operator U, at least for some given time interval At. A negative power is chosen as 
a matter of convenience, recalling that, according to Eq. (6), U~ k = U N ~ k . Fairness of such 
hypothesis shall be justified a posteriori. 

Hypothesis (15) has a direct consequence. If it holds for a given At, then it follows that, 
for any integer n, 

iHAr^ n 



exp 



exp 



h 



iHnAT 

n 



(ir*)' 



u 



-kn 



(16) 



(17) 



and cyclicity of U implies that 



exp 



iHnAT 



h 



jj—kn (modA'') 



(18) 



which obviously means that if Eq. (15) holds for a given time interval, a similar equation 
must hold for all its integer multiples. In view of that we can rewrite our hypothesis as 



g{Ar n ) = U~ k ', k'e[0,N-l], 



(19) 



where 



Ar n = nAr, n — 1,2,3..., 



(20) 



and 



k! = kn (modiV)» 



(21) 



what means that hypothesis (15) refers in fact to an infinite sequence of time intervals and 
not just to a given one. 

A less straightforward consequence can be extracted with the help of the discrete phase 
space representation of the density operator which is, as in the continuous cases, called 
Wigner function and denoted by p w (m,n;t). This function can be immediately calculated 
from 



p w (m, n; t + At) = Tr[G^(m, n)P(t + At)}, 



(22) 



and if one assumes that the hypothesis (15) is true, then (using Eq. (14)) 



p w (m, n; t + At) = Tr[G^{m, n)lJ- k P{t)U k }. 



(23) 



Now, if we substitute the density operator at time t by the expression giving its decompo- 
sition in the operator basis, 



p w (m,n;t + AT) = Tr[&(m,n)U~ k {jj ]T p w (r, s; t)G(r, s) ) U k ] 

r,s=—h 



(24) 



Pw (m, n; t + At) = Tr[ ]T jjV~ l U- j exp 



j,l=-h 



171 



exp 



2ni . ' 
— (mj+nl) 



x 



N 



ITT 



-xz 



exp 



2ni 



r,s=—h x,z=—h 

Using the Weyl commutation relation between V x and U k , 
p w ( mi n;t + AT) = ^ E E E Tr[V- l U-1U'V*\exp 



N 



(xr + zs) 



U k \. (25) 



j,l=—h r,s=—h x,z=—h 



171 



exp 



2ttz 



(mj + nl) 
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cxp 



2ni 



N 



-xk 



p w (r,s;t)exp 



'in 




_N XZ _ 


exp 



2ni 



N 



The trace contribution on the r.h.s. of (26) is N5^5^\ where 5^ is a Kronecker delta 



xr + zs) 



(26) 



modulo N (i.e., it is only different from zero if n = i (mod iV)), and then 



h h 

p w (m,n;t + At) = — ^ ex P 



2ni 



(j(m — r) + /(n — s — k)) 



Pw (r,s;t). (27) 



j,l=—h r,s=—h 

We first observe that the sums over {j, 1} also yield Kronecker deltas modiV, in fact 
N5^ r 8^} k ^ s , and after the sums over {r, s} are performed we finally obtain 



p w (m, n;t + At) = p w (m, n - k;t). 



(28) 



Equation (28) shows that hypothesis (15), for that precise time interval for which it holds, 
leads the system state to a situation equivalent to a shift of k sites along one direction 
in the two-dimensional discrete phase space. Use of equation (19) characterizes a time 
evolution, since for each time interval there will be associated a given shift in phase space. 
A particularly simple and elucidative physical situation where this occurs is presented in 
reference [12]. 

We are yet to determine the physical conditions which allow our hypothesis to be valid. In 
order to achieve that we make use of the discrete phase space representation of the operator 
equation proposed as our hypothesis. The mapped expressions of the individual operators 
read 



(U k )(m, n) = exp 



2ni 



N 



km 



(29) 



(Q(At))(m, n) = g(m, n; At) = exp 



%E m At 



(30) 



respectively, so that on the discrete phase space Eq. (15) will be written as 



exp 



iE m At 

h 



cxp 



2ni 



N 



km 



(31) 



and therefore 



E^= 2 JL km ( mod 27r). (32) 

Equation (32) is in fact a set of N equations (one for each m) which must hold separately 
for a given At and for a given k. It can be still written as 

E m ^ = km + Nf(m). (33) 

where f(m) is an arbitrary integer function of an integer variable. The r.h.s. of Eq. (33) 
is an integer number, whereas the l.h.s is a product of two, in principle, real numbers. It is 
then clear that a necessary but not sufficient condition to Eq. (33) holds is that there must 
be a real number A for which the whole set j^ 1 } is an integer. If we want that condition to 
become sufficient, the set must be additionally a complete set of remainders modulo N [17]. 
Summarizing, Eq. (31) will only have solutions if and only if the spectrum can be written 



as 



E m = hu(km+ Nf(m)), (34) 

where u is an appropriately defined constant. From (32) and (34) it follows that the smallest 
time interval for which G(At) = U~ k is 

^ - m (35) 

and for all integer multiples of At the time evolution operator will be equal to a given power 
of U. Once iV is a prime number, the set of the first N powers of U~ k will be equivalent to the 
set {1, U, U 2 , .., l^ -1 } in some permuted order. Therefore, all powers of U (and consequently 
the respective sites of the discrete phase space) are visited as the time intervals are counted. 

IV. TIME INTERVAL OPERATOR 

We have noticed that, for Hamiltonians which fulfill Eq. (34) and on time intervals 
which are multiples of (35), the time evolution operator is equal to an integer power of the 
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Schwinger unitary operator U. If that is true, the Weyl commutation relation can be written 

as 



exp 



-iH(nAr) 

h 



V j = exp 



V ■' exp 



-iH(nAr) 

h 



(36) 



On the other hand, we can write the Schwinger operator V itself as the exponential of an 
operator if we define 



2tt 



N-l 



(37) 



3=0 



as it can be directly verified that 



exp 



iT(_£/ m _|_ s E m 

h 



v- k \ 



(38) 



once the energy spectrum satisfies (34). Equation (36) therefore can be written as 



exp 



—iHnAr 



h 



exp 



-iTAEj 



exp 



2ni 



■nkjk 



exp 



-iTAEj 



exp 



—iHnAr 



, (39) 



where AEj = E s+ j — E s . 



From Eq. (37) it is not difficult to obtain the eigenstates and eigenvalues of T, namely 



Th)=t,h) = ArZ|t; J >, 



(40) 



and its phase space representative 



(T) (m, n) = Am. 



(41) 



A. Quantum Clock 

Let a given physical system S, described by a Hamiltonian whose spectrum obeys Eq. 
(34), be, at an instant to, in an eigenstate of the time interval operator. Its corresponding 
initial discrete Wigner function is 

p w (m,n;t ) =Tr[G^(m,n)\vi){vi\] (42) 
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which can be easily calculated as 

p w (m,n;to) =6$. (43) 
Using the result of Eq. (28) for the time evolution 

p w (m, n; t + At) = p w (m, n + k; t ) 

it is simple to verify that 

p w {m,n;t + At) = 5 [ ^} +k , 
which can be generalized, through the use of (19) to 

p w (m, n- 1 + jAr) = tfjji-jk- 

This is a quite interesting result. It simply states that, when the system starts its evolution 
at (or visit in a given instant) an eigenstate of the time interval operator, it will be in another 
eigenstate, with no uncertainty, after every time interval At. All eigenstates shall be visited 
in the time interval (iV— 1)At, although, if k ^ 1, they may not appear in the correct order. 
In general the sequence will be 

\vi), \v i+k ), \v i+2 k) |vi+(jv-i)fc), (44) 

where the reader should remind of the modulo N extraction in the subscript of the ket 
labels. In view of that, the next ket in the above list would be 

\vi+Nk) = K), (45) 

and the system would be back to where it started, that is, within such a cycle, the index of 
state ket is a 'counter' of time. 

This strongly motivates us to denote this dynamically determined At (and its first (N— 1) 
integer multiples) as quantum clock time. 



11 



V. CONCLUSIONS 



The basic result that we have achieved here is a direct consequence of the quantum 
mechanics associated to degrees of freedom characterized by a finite dimensional state space 
and from trying to give quantum mechanical meaning to time intervals rather than to time 
itself. We have identified a Hermitian operator, a quantum mechanical observable, which is 
directly related to time intervals. We have denoted this operator as time interval operator, 
and we stress again that it is not directly identified with the continuous usual time parameter 
t. When equations (15) and (34) are simultaneously satisfied we are dealing with physical 
systems such that this operator displays the following properties: 

i) it is a generator of cyclical shifts in energy state space; 

ii) its exponential obeys the Weyl-Schwinger commutation relation together with the time 
evolution operator; 

iii) it is a Hermitian operator whose eigenstates form a set of time intervals with an 
important role in dynamics; 

iv) the eigenstates of the time interval operator are connected with the energy eigenstates 
through a discrete Fourier transform, therefore these two sets have a maximum degree 
of incompatibility. 

The Weyl-Schwinger commutation relation, however, is only fulfilled for some given val- 
ues of time intervals, exactly the eigenvalues of the time interval operator. The time evolution 
for which the Weyl-Schwinger commutation relations hold could then be called stroboscopic, 
in the sense that, in this case, during this time evolution, it always happens that each one 
of the sites of the energy sector of the discrete phase space individually will be visited and 
will exhaust the total probability when the time evolves as multiples of a given time interval 
(which are the eigenvalues of the time interval operator). By its turn, this particular dynam- 
ics must not be confused with that governed by the usual continuous time parameter, which 
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can always be also carried out. Therefore, following the algebraic structure of those opera- 
tors and their corresponding spectra, the equations characterizing the time interval operator 
allowed us to classify some simple physical systems that behave as a quantum clock. 

In this connection, let us discuss some basic operational requirements related to the 
interpretation of the concept of time. To measure time, there must be a reference to a 
nonstationary quantity - or a property - of a physical system, and, in general, this abstract 
quantity admittedly changes continuously with time. From an orthodox quantum mechan- 
ical point of view, a measurement is related to a Hermitian operator, and the process of 
measuring may destroy the clock itself. This consideration illustrates the importance of the 
identification of a Hermitian operator related to time, and why we have called our time in- 
tervals of 'clock times', since a measure would only, in the worst case, 'reset' the clock. We, 
in fact, have not identified a continuous changing property but rather a discrete changing 
property of a physical system to serve as a 'quantum clock'. This comes as a natural feature 
of the finite dimensional spaces we are dealing with. 
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